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Abstract 

A FORTRAN code for generating the leading SU(3) irreducible representation 
(irrep) of ./V identical spin 1 /2 fermions in a harmonic oscillator mean field is in- 
troduced. The basis states are labeled by N - the total number of particles, the 
SU(3) irrep labels (A, fj,), and S - the total spin of the system. The orthonormalized 
basis states have two additional good quantum numbers: e - the eigenvalue of the 
quadruple operator, Qq; and Mj - the eigenvalues of the projection of the total 
angular momentum operator, Jq = Lq + Sq. The approach that is developed can be 
used for a description of nuclei in a proton-neutron representation and is part of a 
larger program aimed at integrating SU(3) symmetry into the best of the currently 
available exact shell- model technologies. 

1 Introduction 

Successful models for describing energy spectra, transition strengths and other nuclear 
phenomena usually work well in one region of a shell but not in others. For example, the 
random phase approximation (RPA) [Q is a reasonable theory for describing properties 
of nuclei near closed shells but fails in mid-shell regions where deformation is the most 
characteristic feature due to the importance of the quadruple-quadrupole (Q ■ Q) inter- 
action [Q] in this domain. For near mid-shell nuclei, algebraic models based on SU(3) [[| 
work best since the basis states are then eigenstates of Q ■ Q. 

Applications of a nuclear shell- model theory always involve three considerations M: 

• Step 1: Selection of a model space based on a simplified Hamiltonian. Frequently 
used approximations include the Simple Harmonic Oscillator Hamiltonian for near 
closed-shell nuclei and a Nilsson Hamiltonian for deformed systems. 

*Paper published in: "Revista Mexicana de Fisica" 44 Suplemento 2, 43-47 (October 1998) 
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• Step 2: Diagonalization of a realistic Hamiltonian in the model space to obtain 
the energy spectrum and eigenstates of the system under consideration. Important 
components of a realistic Hamiltonian include the Q ■ Q and pairing (P) interactions 
as well as single-particle I ■ I and / • s terms. 

• Step 3: Evaluation of electromagnetic transition strengths (E2, Ml, etc.) between 
calculated eigenstates and a comparison with the available experimental data. 

If the model Hamiltonian includes free parameters, for example, single-particle energies 
and/or the strength of the Q ■ Q and P interactions, the procedure is repeated to obtain 
a best overall fit to the experiment data. 

The easiest to use of modern shell-model codes are based on the so-called M-scheme 
logic, namely, model spaces spanned by many-particle configurations (Slater determi- 
nants) with good projection of the total angular momentum |5], 0. Good total angular 
momentum, which is a conserved symmetry due to the isotropy of space, is obtained 
by angular momentum projection. Codes of this type normally yield a good description 
of the single-particle nuclear phenomena; unfortunately, an equally good description of 
collective phenomena within the framework of such a theory is very difficult to achieve. 
On other hand, an SU (3) based shell-model scheme is designed to give a simple interpre- 
tation of collective phenomena. An ideal scenario would incorporate both, allowing the 
Hamiltonian of the system to "choose" which one of the two (or an admixture) is most 
appropriate. 

A project, code named Red Stick for Baton Rouge, is now under development at 
Louisiana State University. Principal authors are Jerry Draayer, Erich Ormand, Calvin 
Johnson, and Vesselin Gueorguiev. The goal of the project is to develop a new M-scheme 
shell-model code that allows the user to include SU (3) symmetry-adapted basis states in 
addition to (the usual) single-particle based configurations. 

In this paper we discuss the first stage of the Red Stick project - a basis generator 
for M-scheme SU(3) shell-model configurations. We begin with a review of SU(3) shell- 
model basics [0, |J. Then we introduce an appropriate single-particle level scheme and 
give matrix elements of the generators of SU(3) in that representation. Next we consider 
the structure of the Highest Weight State (HWS) of an SU(3) irrep, and especially the 
HWS of so-called leading SU(3) irreps. Once a HWS is known, we then go on to show 
how all states of that irrep can be constructed using SU(3) step operators ||. States 
with good projection of the total angular momentum Mj are obtained by considering the 
direct product SU(3) ® SUs(2). Some information about the FORTRAN code that uses 
these techniques will also be presented. 

1.1 5*7(3) Basics 

In this section we review group theoretical concepts that are important to the development 
of the theory and introduce SU (3) conventions adopted in our discussion. We also consider 
the physical reduction SU(3) D 5*0(3) and the canonical reduction SU(3) D U(1)®SU(2) 
with their respective labels. 
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First consider the so-called physical reduction SU(3) D 50(3), which yields a con- 
venient labeling scheme for the generators of SU(3) in terms of SO (3) tensor operators. 
The commutation relations for these SU(3) D SO (3) tensor operators are given in terms 
of ordinary 5*0(3) Clebsch-Gordan coefficients (CGC) j7|: 



Within this reduction scheme, states of an SU (3) irrep (A, fi) have the following labels: 

• (A,/i) - SU(3) irrep labels; 

• I - total orbital angular momentum, which corresponds to the second order Casimir 
operator of 5*0(3); 

• mi - projection of the angular momentum along the laboratory z-axis; 

• k - projection of the angular momentum in a body-fixed frame, which is related to 
multiple occurrences of 50(3) irreps with angular momentum I in the (A,//) irrep. 

Unfortunately, this scheme has only one additive label, namely mi, and in addition, there 
are technical difficulties associated with handling the k label. 

The labeling scheme for this project is the canonical reduction SU(3) D U(l) <g>SU(2) 
where Q is the U(l) generator and the SU(2) generators are proportional to L , Q + 2, 
and Q-2 0. Under the action of the generators of these U(l) and 577(2) groups, the 
remaining four generators of SU(3) transform as two conjugate spin [~] SU{2) tensors 
with e = ±3 values for Q Q . For this scheme, states of a given 577(3) irrep (A, fx) have the 
following labels: 

• (X, fx) - SU(3) irrep labels, 

• e - quadruple moment (Qo), 

• mi - projection of the orbital angular momentum (L ), 

• n p - related to the second order Casimir operator of SU(2), which for symmetric 
(A, 0) irreps is simply the number of oscillator quanta in the (x, y) plane. 

This canonical reduction, SU(3) D U(1)<S> 577(2), has two additive labels, e (Q ) and 
mi (L ) and the allowed values of these labels for fixed SU(3) irrep (A, fi) are given by 0: 




— V2(lm, lm'\lm + m')L m+m i 

— V6(2m, \rd\2m + m')Q m+m t 
3Vl~0(2m, 2m'\lm + m')L m+m i 



(1) 



£ = 2A + / u-3(p + g) 
n p = fi+(p-q) 
mi = rip — 2m 



(2) 



where 0<p<X,0<q<fi, and < m < n p . 
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2 Basis State Generation 



The algorithm we use to generate SU (3) symmetry adapted states with good projection 
of the total angular momentum is described in this section. It consists of four basic 
components: 1) definition of the single-particle levels and matrix elements of the SU(3) 
generators for a given shell; 2) generation of the HWS of SU(3) ® SUs(2) for given 
number of fermions iV and spin S; 3) generation of the other SU(3) states by applying 
step operators to the HWS; and 4) generation of states of good total angular momentum 
projection, Mj. This will be followed by a short discussion of a FORTRAN code that 
uses the algorithm. 

2.1 Single-particle Levels — Ordering Scheme 

Single-particle levels of the rj = 0,1,2,... (s,p, sd, ...) harmonic oscillator shell belong 
to the symmetric (77, 0) irrep of SU(3). Because \i = 0, a typical three-dimensional 
representation of SU(3) basis states, Fig.|Tj, reduces to a special two-dimensional triangular 
shape (e and n p become linearly dependent), Fig.Q. Also, because SU(3) is a compact 
group its irreps are finite dimensional and many-particle (fermion) configurations can 
be conveniently represented as binary strings with a 1 or representing the presence 
or absence of a particle in the corresponding single-particle level. (The latter, together 
with a "sign rule" to accommodate fermion statistics, is a convenient computer native 
implementation of a Slater determinant representation of the basis states.) 

A convenient ordering scheme (which tracks the arrows in Fig.0) is set by requiring a 
simple representation of many-particle configurations with maximum quadruple deforma- 
tion. This objective can be achieved if the states are ordered by e (quadruple moment) 
first and then by mi (projection of the total angular momentum). 

2.2 Action of SU (3) Generators on Single-particle States 

To be able to apply the SU(3) generators on many-particle configurations it suffices to 
know the action of these generators on the single-particle states. The eight generators 
of SU(3) belong to the self-adjoint (1, 1) irrep of SU(3). The operator structure should 
be chosen in the most convenient form for the application under consideration. For the 
present application, this choice is the same as used for the basis states, namely, the 
SU(3) D SU{2) (g> U{1) reduction. The matrix elements of the SU{3) generators can 
be obtained either by using an application of the appropriate Wigner-Eckert theorem or 
by using explicit expressions || for determining the action of the operators on the basis 
states. For computational purposes it is better to adopt a direct solution, one that exploits 
the fact that the action is on a product of single-particle levels each of which belongs to a 
symmetric (77, 0) irrep of SU(3). This allows the matrix elements of the SU(3) generators 
to be calculated using properties of the SU(2) only (Fig||). 

A key feature is the fact that the six non-diagonal generators of SU(3) (recall that 
L and Qq are diagonal) are rising or lowering generators of SU{2) subgroups of SU{3). 
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The three SU (2) subgroups and their respective actions are shown in Fig. |3|. States that 
are collinear with one of the sides of the triangular shape shown on the right in the figure 
form an irrep of the corresponding SU (2) subgroup. 

2.3 Highest Weight States of SU(3)®SU S {2) for Leading Irreps 

So far we have constructed single-particle states and evaluated matrix elements of the 
generators of SU(3) when they act on these states. The next step is to construct many- 
particle HWS of SU(3) ® SUs(2). In the chosen scheme there are seven extreme states 
which correspond to the vertexes of the three-dimensional diagram (Fig.|l]) of a general 
(A, /i) irrep. We are particularly interested in the vertex that has the maximum value for 
the quadruple moment of the system (Fig.[l|). Our HWS is the state with e = 2 A + /x, 
n p = /i, and mi = fx. This HWS (maximum value of mi for maximum e) can be easily 
constructed by ensuring that the action of the SU(3) rising generators annihilate it. 
Indeed, for such a HWS the values of A and fx can be determined from its e and mi labels. 

Selecting the leading (A, fx) irrep (HWS with maximum overall value of e) out of all 
possible irreps of an N fermion system with total system spin S is very simple within the 
chosen scheme. This is because the number of particles with spin up nj and spin down 
rii is uniquely determined by the solution of two linear equations: 

iV = n-f + n± 
2S = nj — nj_. 

Where the second equation expresses the fact that we also require the state to be highest 
weight with respect to SUs(2). Further, maximizing the value of Qq is achieved by filling 
the single-particle states of the (r), 0) irrep (Fig.^J) from bottom to top. The chosen scheme 
ensures that this simple procedure gives maximum values for e and m\. The SU(3) irrep 
labels (A, fx) are obtained by evaluating the quadruple moment Q and projection of the 
angular momentum Lq, as these are additive quantum numbers. 

For example, in the sd-shell, there are six single-particle levels corresponding to the 
(2, 0) irrep of SU(3). The HWS of the leading irrep for N = 6 particles and total system 
spin S = 1 is (3,3) (Fig.fD, whereas for S = the leading SU(3) irrep is (6,0) (Fig.||). 
These many-particle configurations are HWS with respect to SU(3) as well as S77s(2). 

2.4 Generating SU (3) States with Step Operators 

Once we have the HWS of SU(3) ® SUs(2), we can generate any other state of the SU(3) 
irrep (A,/i) by applying step operators similar to those given by Hecht. By using the 
parameterization (|2]) ||, we can identify the corresponding step operators, Fig. |. 

It is important to note that applying p-move or g-move step operators to states on 
the top surface yields other states (or zero) on that same surface. Since the states on 
the top surface are HWS with respect to 577(2) in the SU(3) D SU{2) ® U(l) reduction, 
the m-move step operator is an SU(2) lowering operator which changes the projection of 
the angular momentum, m\. The p-move and g-move step operators can be obtained by 
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imposing the restriction that they generate only transformations within the SU(2) HWS 
space. From an algebraic perspective, the p-move and the m-move operators are linear 
in S77(3) generators while the g-move operator is quadratic. Nevertheless, the state 
generation process can be written in such a way that the g-move operator effectively 
reduces to a linear action. These step operators can also be obtained by a projection 
operator technique JTTJ . 



2.5 Generating States of Good Mj by Lowering the Spin 

Since the action of SU(3) commutes with that of the spin, SUs(2), it is not difficult to 
achieve the final goal of states with good total angular momentum projection, Mj. Recall 
that we started with the HWS of SU(3) ® SUs(2). We then introduced a procedure to 
generate other states of the SU (3) irrep by applying step operators. Each of these states 
remains a HWS of SUs(2). By applying spin lowering operators to SU(3) states with 
mi = Mj — m s , we can generate states with labels: N, S , (A, fx), e, n p , mi, and Mj. 



2.6 About the FORTRAN code 

A FORTRAN code that uses the above algorithm was run on three platforms: Macintosh 
Ilci, Power Macintosh 8500, and SUN workstation. Here we give some of the run-time 
characteristics of the code for shells below the sdg-shelh 

• On a 50 MHz Macintosh Ilci with math co-processor, 3 Mb of memory suffices and 
the longest run time for generating all states of any leading SU(3) irrep was less 
than 5 min. 

• On a 200 MHz Power Macintosh 8500, 4 Mb of memory was sufficient and the 
longest run time for generating all states of any leading SU(3) irrep was less than 
30 sec. 

• On a 200 MHz SUN workstation the longest run time for generating all states of any 
leading SU(3) irrep was less than 15 sec. All states with good Mj for any leading 
irrep of 577(3) ® SUs(2) for the p, sd, and pf shells can be generated in less than 
20 min. 

Beyond the pf -shell (77 > 4), the number of configurations in a single SU(3) state may 
exceed 100,000. The storage requirement must be pushed up accordingly. More storage 
is also needed because beyond the sdg-shell the number of single-particle states exceed 
16 (32 with the spin degree of freedom) which means more than a single word is required 
to represent many-particle configurations. 



3 Summary and Discussion 

An 577(3) basis state generator for producing states with good projection of the total 
angular momentum, Mj, and therefore one that is suitable for integration into modern 
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M-scheme shell-model codes, is suggested as a possible way to obtain simultaneously 
a description of collective and single-particle nuclear phenomena. Indeed, a basis state 
generator that uses the algorithm described above and generates SU(3) states as a sum of 
Slater determinants (many-particle configurations) with good Mj is available. The basis 
state generation procedure consists of four main components as follows: 

• Symmetric irreps of SU (3) are chosen as single-particle shell-model states, 

• many-particle state is constructed as HWS of SU(3) ® SUs(2), 

• SU(3) step operators are applied to obtain other states within an SU(3) irrep, 

• and finally, states with good Mj are then obtained by SUs(2) lowering operators. 

There are additional things that need to be done, especially the generation of other 
then leading irreps of SU(3). This can be achieved by employing an algorithm for gen- 
erating non-leading HWS that is part of an SU(3) reduced matrix element package 
or perhaps by using another more direct approach tailored to the present application. 
The generation of other HWS is important if one wishes to include states that are not 
maximally deformed in their intrinsic configuration in a calculation. For example, these 
will be important if the non-Q ■ Q parts of the interaction play a significant role. 

The generation of good proton-neutron SU(3) and Mj states can be achieved in a 
variety of ways. One scheme, the so-called strong coupling limit uses standard SU(3) 
coupling procedures to couple the proton and neutron irreps to final irreps: proton 
states, \N n S" 1 (A 71 ", / u 7r ) e n , n*, mf, MJ >, are coupled with neutron states, \N" S u (A", /i v ) 
e u ,np,mf , Mj >, to final states \N S (A,/i) e,n p ,mi, Mj >. Another procedure would 
be to extend what was done here for SUs(2) to SU(4) reduced with respect to SUs(2) <g> 
SUt(2), the supermultiplet scheme. The latter is not simple because the SU(4:) HWS is 
not the one of primary interest in nuclear physics and therefore it must be reduced, which 
requires procedures like those developed here for a non-canonical reduction of SU(A). 

The underlying goal is to be able to carry out shell-model calculations using M-scheme 
techniques with SU(3) symmetry adapted basis states included in the basis. The SU(3) 
basis state generator described here provides good start in that direction, and therefore 
the beginning of an symmetry adapted M-scheme code that is an integral part of our Red 
Stick shell-model project. 
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Figure 2: Ordering scheme of the single-particle levels. 



9 



Figure 3: The diagram on the right shows the action of the six non-diagonal generators 
of SU(3). The vertical solid lines represent the action of the SU{2) subgroup that enters 
the SU(3) D SU(2) ® U(l) chain. The diagram on the left shows that applying an SU(3) 
generator to a single-particle state results in another single-particle state. 



e = -2 m,= -2 

e = -2 m, = 

e = -2 t i > m,= 2 

e= 1 4> m,= -l 

e = 1 4> ^ 111,= 1 

e = 4 ^ m, = 

m, = 3 = (j, and e = 9 => X = 3 

Binary Form: 0000...0 001111 000011 



Figure 4: Highest weight state of the leading (3, 3) irrep for N = 6 and S — 1 in the 
sd-shell. 
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E = -2 ■>>, = -2 

e = -2 in, = 

£ = -2 m,= 2 

£ = 1 ' 111, = -1 

E = 1 4= m,= 1 

E = 4 't' in, = () 

m, = O = (X and £ — 12 => ^ = 6 
Binary Form: 0000...0 OOOlll OOOlll 

Figure 5: Highest weight state of the leading (6, 0) irrep for N = 6 and S = in the 
sd-shell. 



\X q-moves 




Figure 6: Action of the SU(3) step operators. 
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